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ABSTRACT 

We report the discovery of multiple condensations in the prestellar core can- 
didate SMM 1A in the R CrA cloud, which may represent the earliest phase of 
core fragmentation observed thus far. The separation between the condensations 
is between 1000 and 2100 AU, and their masses range from about 0.1 to 0.2 M Q . 
We find that the three condensations have extremely low bolometric luminosities 
(< 0.1 Lq) and temperatures (< 20 K), indicating that these are young sources 
that have yet to form protostars. We suggest that these sources were formed 
through the fragmentation of an elongated prestellar core. Our results, in con- 
cert with other observed protostellar binary systems with separations in the scale 
of 1000 AU, support the scenario that prompt fragmentation in the isothermal 
collapse phase is an efficient mechanism for wide binary star formation, while the 
fragmentation in the subsequent adiabatic phase may be an additional mecha- 
nism for close (< 100 AU) binary star formation. 

Subject headings: binaries: general — ISM: clouds — ISM: individual (R Corona 
Australis, R CrA SMM 1A) — stars: formation 



1. INTRODUCTION 

The origin of binary stars is still a puzzle in our understanding of star formation. Nu- 
merous theoretical simulations support the hypothesis that the fragmentation of collapsing 
molecular cloud cores, promoted by either rotation or turbulence, is the main mechanism for 
the formation of binary/multiple stellar systems (see reviews by Bodenheimer et al. 2000, 
Tohline 2002, and Goodwin et al. 2007). Nevertheless, many key questions concerning this 
fragmentation process, e.g., when exactly does the fragmentation start, are still under debate 
(see Tohline 2002). Although it is generally assumed that cloud cores do not fragment during 
the free-fall collapse phase, several groups (e.g., Bate & Burkert 1997; Tsuribe & Inutsuka 
1999; Boss et al. 2000) found that fragmentation can occur at the end of the isothermal 
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phase (generally referred to as prompt fragmentation), while others (e.g., Truelove et al. 
1997, 1998; Machida et al. 2005, 2008) argue that the isothermal gas is stable against frag- 
mentation and found that fragmentation only occurs during/after the adiabatic phase (see 
Figure 9 in Andre et al. 2009 for a discussion of the different evolutionary phases of core 
collapse) . 

On the observational side, a handful of young protostellar (i.e., Class 0) binary sys- 
tems have been found (e.g., Looney et al. 2000; Launhardt 2004), and there are increasing 
kinematics studies of binarity in the protostellar phase (e.g., Chen et al. 2007, 2008, 2009; 
Volgenau et al. 2006; Maury et al. 2010). However, the number of observed and well-studied 
protostellar binary (protobinary) systems is still small, and these systems have already gone 
through the fragmentation process. The observational link between the initial conditions in 
a fragmenting cloud core and the final stellar systems formed therein is still missing. It is 
therefore critical to observe, at high angular resolution, more dense cores in nearby molecular 
clouds to search for the earliest phase of fragmentation and study in detail their properties, 
in order to put direct constraints on fragmentation models. 

In this Letter, we present Submillimeter Arra}0 (SMA; Ho et al. 2004) dust continuum 
observations toward the R Corona Australis (RCrA) region. At a distance of ~170pc 
(Knude & H0g 1998), the CrA dark cloud is one of the nearest star- forming regions. As the 
most active star formation site in this cloud, the R CrA region has been extensively observed 
in the past two decades (see review by Neuhauser & Forbrich 2008, and references therein). 
Using SCUBA on the James Clerk Maxwell Telescope (JCMT), Nutter et al. (2005) found 
a prestellar core candidate termed SMM 1A in the R CrA region, which is elongated in the 
east-west direction and has a gas mass of ~1OM , and an effective radius of ~3000AU. 
With a maximum velocity dispersion of about 0. Skins" 1 (Harju et al. 1993), the SMM 1A 
core is gravitationally bound. Infall motions in this region of the cloud further confirm that 
this is a prestellar core (Groppi et al. 2004). In this Letter, we report the discovery of 
a multiple system within the SMM 1A core, based on high angular resolution millimeter 
observations. This may represent the earliest phase of core fragmentation observed thus far. 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. OBSERVATIONS AND DATA REDUCTION 

The R CrA region was observed with the SMA on 2006 August 20 in the compact con- 
figuration. Six antennas were used in the array, providing baselines from 5kA to 52 kA at 
220 GHz. The SMA primary beam is about 55" at this frequency. Two positions next to each 
other at (R.A., decl.)j2ooo =( 19:01:53.3, -36:57:21.0) and (R.A., decl.) J2 ooo =( 19:01:56.4, 
—36:57:27.0) were observed. The digital correlator was set up to cover the frequency rang- 
ing 216.6— 218.8 GHz and 226.6— 228.8 GHz in the lower and upper sidebands, respectively. 
The 1.3 mm dust continuum emission was recorded with a total bandwidth of ~ 3.3 GHz 
(~ 1.8 GHz USB and ~ 1.5 GHz LSB). System temperatures for RCrA ranged from 100 to 
280 K (depending on elevation), with a typical value of ~200K. 

The visibility data were calibrated with the MIR package (Qi 2005), with quasars 3c454.3 
and 1924-292 as the bandpass and gain calibrators, respectively. Uranus was used for ab- 
solute flux calibration, from which we estimate a flux accuracy of ~20%, by comparing 
the measured quasar fluxes with those listed in the SMA calibration database. The cali- 
brated visibility data were imaged using the Miriad toolbox (Sault et al. 1995). The SMA 
synthesized beam size at 1.3 mm dust continuum, with robust uv weighting 0, is 5.7" x 2.3". 

3. RESULTS 

Figure la shows the SMA 1.3mm dust continuum image of RCrA, overlaid with the 
SCUBA 850 /im dust continuum contours (from Groppi et al. 2007). The northern part of 
this image shows clear detection of the dust continuum emission associated with the Herbig- 
Ae star RCrA and protostellar cores SMM2, SMM IB and SMM1C — these last two are 
also referred to as SMAl and SMA2 in Groppi et al. (2007; see also Choi et al. 2008). To 
the south, the SCUBA 850 /iin image shows the SMM 1A core, which is elongated in the 
east- west direction (see Fig. la). At higher angular resolution, the SCUBA 450 /im image 
in Figure lb shows that the SMM 1A core is clearly resolved into two sub-cores (see also 
van den Ancker 1999). At even higher resolution, the SMA 1.3 mm continuum observations 
reveal that the SMM 1A core is divided into three condensations, aligning roughly from east 
to west, which we refer to here as SMMlA-a, SMM lA-b, and SMM lA-c (Fig. lb). All are 
detected with signal-to-noise ratio of 4 or more. Source SMM lA-a coincides with the eastern 
sub-core observed in the SCUBA 450 jim image, while SMM lA-b and SMM lA-c are coin- 
cident with the western sub-core (Fig. lb). The single-dish and interferometer observations 
are in general agreement with each other, which indicates that the three continuum sources 
detected with the SMA are neither artifacts nor are due to noise in the interferometer im- 
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agc|. Also note that the SiO (5-4) line and several other molecular lines, e.g., DCN(3-2), 
were included in this SMA correlator setup, but no line emission was detected from the three 
condensations. 

Figure 2 shows the mid-infrared images of the R CrA region, taken by the Spitzer Space 
Telescope. The MIPS 1 band image at 24 /im (Fig. 2a) clearly shows infrared emission at 
the position of cores SMM IB (associated with the infrared source IRS 7B, see Wilking et 
al. 1997), SMM 2, and the Herbig stars RCrA and TCrA. There is also 24 /im emission 
from the SMM1C core, which may be associated with the radio source B9 centered on the 
core (see Choi et al. 2008), and the infrared source IRS 7A located ~5" to the south. The 
MIPS 2 band image at 70 /im (Fig. 2b) has a lower angular resolution and is dominated by 
the emission from the SMM 1C/IRS 7A region. The MIPS 3 (160 /xm) image of RCrA is 
extremely saturated, and is thus not shown here. None of the Spitzer images at bands from 
3.6 to 70 /im shows compact infrared emission in the SMM 1A core, which suggests that the 
core is extremely cold. 

From the SMA 1.3 mm images, we derived the positions, fluxes, and (deconvolved) 
FWHM sizes of the three sources in R CrA SMM 1A using a mult i- component Gaussian 
fitting routine (see Table 1). Angular separations between sources SMM lA-a and -b and 
SMM lA-b and -c are 12 / .'2±0'/3 and 5'.'8±0'.'3, corresponding to projected separations of 
~2100AU and ~ 1000 AU (at a distance of 170 pc), respectively. Assuming that the 1.3 mm 
dust continuum emission is optically thin, the total gas mass (M gas ) of the three sources 
is calculated with the same method as described in Launhardt & Henning (1997). In the 
calculations, we adopt a dust opacity of Kd = 0.5 cm 2 g -1 (Ossenkopf & Henning 1994), a 
typical value for dense and cold molecular cloud cores, and a mass-averaged dust temperature 
of ~ 18 K (see below). The gas masses of the three sources, derived from the SMA 1.3 mm 
dust continuum observations, range from ~0.10 to ~O.23M (see Table 1). With these 
mass and the FWHM size of the condensations, we estimate their average densities to be 
in the order of 10 7 to 10 s cm -3 . Note that these are average densities and the local peak 
densities in these sources may be orders of magnitude higher. These high densities, as well 
as the infall motions detected in this region (see Groppi et al. 2004), suggest that the three 
condensations in the SMM 1A are on their way to form low-mass stars. Considering that the 
SMM 1A core has a mass of approximately 10 M , and assuming a core-to-star efficiency of 
~ 30% (Evans et al. 2009), it is probable that a multiple stellar system with a total mass of 



2 The three sources in SMM 1A were not detected in the SMA 1.1 mm observations in Groppi et al. (2007), 
probably due to the reasons: (1) the three sources lie near the edge of their primary beam (~ 46" at 270 GHz) ; 
(2) more flux was resolved out in their extended configuration observations; and (3) the 1.1 mm observations 
lack the sensitivity to detect these weak sources. 
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3 M Q will eventually form in SMM 1A. 

4. DISCUSSION 

4.1. The Earliest Phase of Binary/Multiple Star Formation 

Table 1 lists the 450 pm and 850 fim fluxes estimated from the SCUBA images, and the 
800 pm fluxes adopted from van den Ancker (1999). Since sources SMM lA-b and SMM 1A- 
c are not resolved in the JCMT observations, we treat them as one sub-core here. The 
SCUBA 450 ^m fluxes are derived by fitting two-dimensional Gaussians toward the two sub- 
cores (SMM lA-a and SMM lA-b+c), while the 850 fluxes are estimated from the fluxes 
enclosed within roughly one beam around each sub-core. Based on the SMA 1.3 mm and 
JCMT submm data points, as well as the 3 a upper limits in the Spitzer MIPS images, we 
construct the spectral energy distributions (SEDs) of the two sub-cores (see Figure 3). In 
order to derive luminosities and temperatures, we first interpolate and then integrate the 
SEDs, always assuming spherical symmetry. Interpolation between the flux densities is done 
by a x 2 single-temperature grey-body fitting to all points (including the upper limits). The 
resulting bolometric luminosities are <O.O8L for SMM lA-a and <O.O9L for SMM lA- 
b+c. The dust and bolometric temperatures are between 17 and 19 K for the sources. 

The low luminosities and temperatures, as well as the fact that no compact infrared 
emission is detected from these sources in the Spitzer images, resemble the typical properties 
of prestellar cores (see review by Andre et al. 2009). In the VLA survey of the RCrA region, 
no radio source is detected at the positions of the three condensations presented here (Choi 
et al. 2008). In contrast, the protostars SMM IB (Class0/I) and SMM 1C (ClassO) to the 
north of SMM 1A are bright and centrally peaked at millimeter wavelengths (Fig. la), and are 
associated with hard X-ray (Forbrich & Preibisch 2007), infrared (Fig. 2a), and centimeter 
radio counterparts (Choi et al. 2008). Therefore, although the chemical and kinematic 
properties of SMM 1A are still poorly known, the observations thus far suggest that the 
three sources in SMM 1A are in an earlier evolutionary stage than that of Class protostars. 
As predicted by theoretical studies, there are different phases in the evolution from an 
initial isothermal (10 K) prestellar core to a Class protostar (Andre et al. 2009). The low 
luminosities (< 0.1 L ) and temperatures (close to 20 K) of the three sources in SMM 1A are 
similar to the predicted properties of so-called first hydrostatic (adiabatic) cores, which are 
formed after the isothermal collapse of prestellar cores (e.g., Masunaga et al. 1998; Chen et 
al. 2010). However, it must be noted that large uncertainties remain in our estimates due to 
the limited observations available. Further high angular resolution observations at different 
wavelengths are needed to constrain the SEDs of the three sources in order to address more 
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precisely their evolutionary statuses. Regardless of this, it is fair to say that the R CrA 
SMM1A multiple system represents the earliest phase of core fragmentation observed thus 
far. 

4.2. Prompt (Isothermal) Fragmentation vs. Adiabatic Fragmentation 

Prompt fragmentation of rotationally flattened cloud cores with initially flat density 
profiles, immediately after a phase of free-fall collapse, has long been suggested as one of the 
most efficient mechanisms for binary star formation (Tohline 2002). Since this fragmenta- 
tion is expected to occur at the end of the isothermal collapse phase, it is also referred to 
as isothermal fragmentation. The separation of the fragments formed via prompt fragmen- 
tation has a scale of 10 2 — 10 4 AU, which corresponds to the Jeans length in this phase. In 
contrast, other groups found that isothermal gas is stable against fragmentation and that 
core fragmentation can only occur during (and after) the adiabatic phase (e.g., Truelove et 
al. 1998; Machida et al. 2008). Because this fragmentation occurs after the formation of the 
first adiabatic core, it is also called the first core fragmentation (Machida et al. 2005). Here, 
the separation of the fragments formed generally has a scale of 3— 300 AU (Machida et al. 
2008). Clearly, observing the earliest phase of core fragmentation is the most effective way to 
discriminate between prompt (isothermal) and adiabatic fragmentation scenarios. If prompt 
fragmentation does occur in the isothermal phase, we may expect to find a binary/multiple 
system consisting of adiabatic cores with the separation of ~ 1000 AU. 

As shown in the SCUBA images (see Figure lb), the RCrA SMM1A core is highly 
elongated with an aspect ration of roughly 3, showing a morphology similar to the bars or 
filaments seen in the simulations of prompt (isothermal) fragmentation (see, e.g., Bate & 
Burkert 1997; Tsuribe & Inutsuka 1999; Boss et al. 2000). From the mass and radius of 
the SMM1A core (see Nutter et al. 2005), the average volume density is calculated to be 
1—2 x 10 7 cm -3 . The corresponding Jeans length is ~ 1000 AU for a temperature of ~ 10 K. 
The separations among the three sources found in the SMM 1A core are also consistent with 
this estimated Jeans length. Moreover, our results from the SED fitting (see above) suggest 
that the three sources are at an evolutionary stage earlier than that of Class protostars, 
perhaps in the first (adiabatic) core stage. All these results appear to support the view that 
the three sources are formed through the prompt fragmentation of an elongated collapsing 
core during the isothermal phase. In this scenario, we would expect the rotational axis of 
SMM 1A to be perpendicular to the direction of the core's major axis, which would be roughly 
perpendicular to the large scale rotation observed in the R CrA region by Groppi et al. (2004). 
This is not necessarily inconsistent with our picture as other cores show drastic differences 
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in the rotation axes measured at different scales (e.g., Ohashi et al. 1997). We also note that 
so far most observed protostellar binary/multiple systems (e.g., L723, Launhardt 2004; or 
NGC2264 D-MM1, Teixeira et al. 2007) have separations in the scale of 1000 AU, which also 
supports the scenario of prompt fragmentation, although these protostellar systems have had 
at least a few 10 4 yrs to evolve from their prestellar stage and thus their separations may 
not represent real initial conditions. 

On the other hand, as introduced above, prompt (isothermal) fragmentation has trou- 
ble accounting for the close binary systems with separations < 100 AU. The fragmentation of 
cores at the adiabatic phase is still an attractive mechanism for close binary star formation. 
Considering that the binary separation distribution in (pre-) main sequence stars peaks at 
~30AU, we actually would expect to find frequently close protobinary systems with sepa- 
rations of < 100 AU to support the adiabatic fragmentation scenario. In practice, however, 
the number of this kind of systems is very small. This observational result can be under- 
stood by the fact that current large millimeter interferometers (e.g., SMA and IRAM-PdBI) 
normally reach 1-2" angular resolution under general conditions, and thus mostly resolve 
wide protobinary systems with separations of 200 AU or more in nearby star-forming clouds. 
We believe that routine observations at angular resolution better than 0.5" will reveal more 
close protobinary systems, which will then allow us to study in detail a statistically signifi- 
cant number of protostellar binary/multipe systems with a wide range of separations (from 
-10AU to 10 4 AU). 



5. SUMMARY 

We present SMA 1.3 mm dust continuum observations toward the RCrA region. The 
1.3 mm dust continuum emission is detected from dense cores SMM IB, SMM 1C, and SMM 2, 
and the Herbig-Ae star R CrA in this region. We discover within the prestellar core candidate 
SMM 1A a multiple system made up of three condensations with masses in the range of 0.1 
to 0.2 M . The angular separations between the three new sources are ~12.2" and ~5.8", 
corresponding to the projected separations of ~ 2100 AU and ~ 1000 AU, respectively. Based 
on SMA observations and complementary data from the JCMT and Spitzer telescopes, we 
construct their spectral energy distributions (SEDs) and find that all three sources have 
extremely low bolometric luminosities (<0.1L Q ) and temperatures (<20K). We suggest 
that the three condensations in SMM 1A resulted from the fragmentation of an elongated 
core in the isothermal phase. The RCrA SMM 1A system appear to be the earliest phase of 
low-mass core fragmentation observed thus far. 
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Table 1. SMA 1.3 mm and JCMT submm dust continuum results of RCrA SMM 1A. 



Source 


R.A. a 


Dec." 


S„(1.3mm)» 


FWHM" 


Mgas 6 


S„(450 /an) 


S„(800 (im) 


S„(850 urn) 




(J2000) 


(J2000) 


[mjy] 


maj.xmin. 


[M e ] 


[Jy] 


[Jy] 


[Jy] 


SMM lA-a 


19:01:55.86 


-36:57:48.8 


46±9 


5'/ 6x5'.' 2 


0.10±0.02 


5.0±2.0 


1.5 


1.0±0.2 


SMM lA-b 


19:01:54.87 


-36:57:45.2 


64±13 


6'.' 8x6'.' 3 


0.13±0.02 


4.6±1.8 C 


0.8 C 


0.7±0.2 C 


SMM lA-c 


19:01:54.56 


-36:57:40.6 


110±20 


7'.'8x3'.'l 


0.23±0.05 


4.6±1.8 C 


0.8 C 


0.7±0.2 C 



a Ccntcr position, flux, and FWHM size of the SMA 1.3 mm dust continuum sources derived from the Gaussian fitting. 

b Total gas mass derived from the SMA 1.3mm dust continuum observations (see text for the dust opacity and temperature used). 

c Combincd fluxes for both sources SMMlA-b and SMMlA-c. 
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Fig. 1. — (a) SMA 1.3 mm continuum data (solid contours) overlaid on the JCMT/SCUBA 850 /im con- 
tinuum map (dashed contours) of the RCrA region. The SMA data contours correspond to 3, 5, 9, 14, 
20cr levels and then increase in steps of 10a (1 a ~ 6.2 mJy beam -1 ). The SCUBA 850 yum contours begin 
at 30% and then increase in steps of 20% of the peak emission (~ 2.4 Jy beam -1 ). The synthesized SMA 
beam (5.7" x 2.3") and the HPBW of the SCUBA (~ 14") arc shown as grey ovals in the bottom right 
corner, (b) SMA 1.3 mm dust continuum image (solid contours) overlaid on the JCMT/SCUBA 450 ^im 
continuum map (dashed contours). Here, the SMA data start at the 3a level and increase in steps of la. 
The SCUBA 450 /im contours (dashed lines) begin at 75% and then increase by a step of 5% of the peak 
emission (~ 10.2 Jy beam -1 ). The synthesized SMA beam and the HPBW of the SCUBA map (-8") are 
shown in the bottom right corner. 
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Fig. 2. — (a) Spitzer MIPS 1 (24 fim) image of the R CrA region, with contours showing 
the SCUBA 450 /imcontinuum emission. Contours levels begin at 20% and increase by steps 
of 10% of the peak emission (~ 10.2 Jybeam -1 ). (b) Spitzer MIPS 2 (70 fim) image of the 
R CrA region, with contours showing the SCUBA 850 ^m continuum emission. Contour 
levels begin at 20% and increase in steps of 10% of the peak emission (~2.4 Jybeam -1 ). 
In both images the crosses mark the positions of the SMA 1.3 mm dust continuum sources, 
while star symbols show the positions of the Herbig stars RCrA and T CrA, and the infrared 
source IRS 7A (the MIPS 1 image is saturated at the position of these three sources). 




Fig. 3. — Spectral energy distributions of the RCrA SMM lA-a and the SMM lA-b+c sub- 
cores. The Spitzer 70 ^m and 160 /xm points, and the SMA 1.3 mm point represent upper 
and lower limits, respectively (indicated by the arrow). Solid lines show the best-fit for 
all points at A > 70 jum using a greybody model S v = B V (T A )(1 — e~ Tv )£l, where B U (T A ) is 
the Planck function at frequency v and dust temperature Td, r v is the dust optical depth 
as a function of frequency, and f2 is the solid angle of the source. The fitting results are 
summarized in the plot. 



